Electron self-trapping in intermediate-valent SmBg 
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SmBg exhibits intermediate valence in the ground state and unusual behaviour at low temperatures. 
The resistivity and the Hall effect cannot be explained either by conventional s/-hybridization or 
by hopping transport in an impurity band. At least three different energy scales determine three 
temperature regimes of electron transport in this system. We consider the ground state properties, 
the soft valence fluctuations and the spectrum of band carriers in n-doped SmBe. The behaviour of 
excess conduction electrons in the presence of soft valence fluctuations and the origin of the three 
energy scales in the spectrum of elementary excitations is discussed. The carriers which determine 
the low-temperature transport in this system are self-trapped electron-polaron complexes rather 
than simply electrons in an impurity band. The mechanism of electron trapping is the interaction 
with soft valence fluctuations. 
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I. INTRODUCTION 



Samarium hexaboride is the first compound in which 
the phenomenon of intermediate valence (IV) has been 
seen directly by X-ray absorption [Q, but the theory 
which explains all properties of this compound within a 
unified physical model remains elusive in spite of unceas- 
ing efforts of theoreticians and experimentalists during 
the recent thirty- five years. 

The conventional description of the electronic struc- 
ture of SmBg in terms of a two-band promotion model 
results in a semiconductor- like spectrum with a gap A. 
This gap appears as a result of on-site hybridization be- 
tween the narrow band formed by electrons from the 
samarium 4/-shell and the wide conduction band formed 
by boron p-states and samarium 5(i-states. The theory 
explains, at least phenomenologically, the main body of 
available experimental data However, the most in- 
triguing properties of this material (as well as other IV 
semiconductors from the same family, i.e. SmS, TmSe, 
YbBi2) cannot be described adequately within a frame- 
work of such a crude phenomenological picture. Among 
these properties are the very nature of the IV ground 
state, the origin of slow valence fluctuations, and the 
low-temperature anomalies of various physicalproper- 
ties. Recent experimental studies of transport y-g and 
optical [^j-^ properties revealed the existence of several 
energy scales in the low-energy excitation spectra and 
several temperature regimes in the low-temperature elec- 
tron kinetics. 

It is found that at least three different activation en- 
ergies influence the behaviour of electrons. The biggest 



energy gap A 



opt 



14 — 20 meV is observed in the fre- 



quency dependent conductivity and dielectric permittiv- 
ity. Apparently, just this value should be ascribed to the 
hybridization gap A of the two-band model mentioned 
above. A substantially smaller value of Aact ~ 3 — 5 meV 
is seen in the low-temperature resistivity and Hall effect 
measurements which display an activation type temper- 
ature dependence in the interval 15 > T > 6 K. This gap 
was also measured in tunnelling experiments Below 
6K the resistivity is nearly temperature independent, and 



such behavior seems to be consistent with hopping trans- 
port in an impurity band separated by the gap Aact from 
the bottom of conduction band. However, there is still 
no experimental consensus about the energy scales which 
control hopping. In ref. H] part of the low-T interval is 
presumed to be described by Mott's T^^^ law for variable- 
range hopping with a scaling energy of Tq = 53 K |^ . In 
another series of samples an activated term in the con- 
ductivity was observed with an activation temperature 
Ta = 2.68 K for T < 3i4: @, although the derivation of 
an activation energy of the order of several K from mea- 
surements made in the same temperature interval looks 
not too trustful. Besides, the pressure dependence of the 
residual resistivity Q] is extremely strong, and it hardly 
can be fitted in a picture of non-interacting electrons in 
an impurity band. In any case, there is no room for ad- 
ditional excitation branches in the mean-field two-band 
theory, and adding extra localized states in the gap docs 
not improve the situation. The standard hybridization 
model seems to be too simplistic, and we believe that 
the generic properties of the intermediate valence state 
can be described only within a framework that goes be- 
yond the mean field approximation. 

In the present paper we offer a description of the low- 
energy spectrum of IV SmBg by treating the phenomenon 
of intermediate valence in rare-earth semiconductors in 
terms of an excitonic dielectric state. This approach was 
offered nearly two decades ago and its effectiveness 
was demonstrated later in explanation of anomalies in 
the vibration spectra of IV semiconductors jl^. Fea- 
tures of excitonic instability were seen also in studies of 
the dielectric-metal transition under pressure in the com- 
pounds TmSei_a;Sa; [Q. It is interesting that recently 
the dramatic observations of a ferromagnetic phase in 
n-doped divalent hexaborides CaBg and SrBg [13| were 
interpreted in terms of conventional excitonic instability 
of a Keldysh-Kopayev- Volkov-Rusinov type jl^] . 

The theory of intermediate- valence excitonic dielectric 
emphasizes the role of soft singlet excitonic states in the 
formation of the ground IV state and valence fluctuations 
in the excitation spectrum. The IV state is described as 
a mixture of singlet ''Fq states of divalent Sm(f^) and 
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the bound electron-hole pairs f^b where 6 is the state 
of an electron promoted from the /-shell to p-orbitals 
spread over neighbouring boron sites but having the same 
symmetry as the f-electron in a central site. In some sense 
these states are electron-hole analogs of Zhang-Rice two- 
hole states known in the theory of low-energy states of 
Cu02 planes in high-Tc perovskites [ p^ . 

Here we consider the case of doped n-SniBg material. 
We study the behaviour of excess conduction electrons in 
the presence of soft valence fluctuations and discuss the 
origin of three energy scales seen in different tempera- 
ture regimes. We show that the carriers which determine 
the low-temperature transport properties in this system 
are self-trapped electron-polaron complexes rather than 
simply electrons in an impurity band. The mechanism 
of electron trapping is the interaction with soft valence 
fluctuations. 



II. GROUND STATE OF INTERMEDIATE 
VALENCE SEMICONDUCTOR 

SmBe together with "golden" SmS is considered to be 
a classical example of a non-magnetic intermediate va- 
lence (IV) semiconductor (see j|] and references therein) . 
Both of these compounds possess a singlet ground state 
with intermediate valence v +2.6 . The main difference 
in the properties of these two compounds is that SmS 
transforms into an IV semiconductor only at finite pres- 
sure (or under the chemical pressure of rare earth ions of 
smaller radius substituting for Sm) whereas SmBg pos- 
sesses intermediate valence and concomitant anomalous 
properties at ambient pressure. 

Apparently, differences in the electron band spectra of 
these two materials are the eventual source of differences 
in their properties. SmS at ambient pressure is known 
to be a normal semiconductor with a relatively wide gap 
in the energy spectrum A = 0.23eV. This gap divides 
the filled /-levels of Sm^+(/®) ions from the bottom of 
conduction band formed mainly by Sm d-states. The 
4/-levels, in turn, form a nearly dispersionless valence 
band within a wide forbidden gap Eg ^ A between the 
conduction band and conventional valence band formed 
mainly by p-electrons of the chalcogen sublattice. At fi- 
nite pressure the gap A is suppressed, but instead of col- 
lapsing to zero it transforms into a microgap A^ at some 
critical pressure Pc, and the material acquires the prop- 
erties of an IV semiconductor. On the other hand, the 
"parent" divalent hexaborides without a 4/-shell (CaBg, 
SrBg), as well as trivalent LaBg with an empty 4/-shell, 
have a band structure without a gap in the actual energy 
interval. The density of states near the Fermi level in hex- 
aborides is predetermined by a single degenerate band 
with a minimum in the X-point of the Brillouin zone. 
This band is formed mainly by p-electrons of the boron 
sublattice with a small admixture of d-states of the cation 
sublattice. It is nearly empty in divalent semimetallic 
hexaborides |l6|jr^ ] and more than half-filled in trivalent 
LaBg Thus, one can expect that the /-level which is 
responsible for the intermediate valence of SmBg should 



cross this simple band. Indeed, conventional band cal- 
culations [p^ give a band structure compatible with this 
presumption. 

Despite the differences between the parent (unhy- 
bridized) spectra of SmS and SmBg, the decisive simi- 
larity of these two cases is the closeness or overlapping 
of the highly correlated 4/-levels of Sm^+(4f^) and the 
nondegcnerate conduction band. In SmBg these states 
overlap at ambient pressure, whereas in SmS the close- 
ness is achieved under external pressure which induces 
the excitonic instability |]l^ at Pc- The closeness re- 
sults in promotion of electrons to delocalized band states. 
Crucial is the fact that the number y of promoted elec- 
trons is the deviation from the divalent state in both 
Sm(2+y)+S and Sm^'^+y^+Be IV semiconductors. Accord- 
ing to the conventional band scheme which treats the //- 
interaction in a self-consistent mean-field approximation, 
the /p-mixing results in the appearance of a hybridiza- 
tion gap A in case of SmBg and a pseudogap in the mixed 
valent phase of SmS. To take into account strong on-site 
correlations, this scheme should be modified. 

The starting Hamiltonian for treating IV semiconduc- 
tors is the Anderson lattice Hamiltonian supplemented 
by an interband Coulomb interaction (see, e.g., pO|): 

iJ = ^ EkaL^k^ + XI -EA|mA)(mA| 

kcr m A=0,{r/j} 

+ [^mka,rA.aL|mrM)(mO| + i7.c. 

+ X C^mkk'nm/ako.ak'cr- (1) 
jkk'cr 

Here ek is a simple spin-degenerate dispersion law for 
the band electrons (we assume for the sake of simplic- 
ity that these electrons come from the boron 2p-shells). 
The states |mA) of the Sm ion at site m are represented 
by two configurations, |mA) = |mO) = |m''_fo) for the 
divalent state 4f^ and |mA) = \mTfi) for the trivalent 
state 4f^. r stands for the muhiplet ^-^5/2 of the 4f^ 
configuration which splits into a Fy/i doublet and a FgyU 
quartet by the cubic crystal field Acf', M enumerates 
the states of corresponding irreducible representations, 
and rijnf = ^ ^ \mr fi) {mT fi\ is the /-electron occupa- 
tion number. The hybridization interaction describes the 
promotion of electrons from the /-shell to the conduc- 
tion band accompanied by a change of atomic config- 
uration. The strong interband Coulomb interaction is 
given by C/^kk' ■ T^^^ hybridization matrix element is 
VmkaXt^ = (kcr,mF/x|i?|mO). At T < Acf only the 
lowest doublet state F7 should be taken into account, so 
one should deal with hybridization of two doubly degen- 
erate states. 

The strong intrashell interaction is inserted in the con- 
figuration change operators — |mA)(mA'|. To 
make the hybridization problem solvable the nondiag- 
onal operators are usually represented in a factorized 
form a:™ = fmp.bln (see, e.g., where /m^ and bl^ 
describe auxiliary fermion and boson fields which corre- 
spond to charge and spin degrees of freedom. This proce- 
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dure gives reasonable results in the case of nearly integer 
valence, but in the IV state the spin charge separation 
procedure can hardly be useful even as a starting approx- 
imation. We prefer to use another approach |22|] which 
seems to be adequate in the nonmetallic case when one 
can hope that hybridization will result in the separation 
of strongly correlated and weakly correlated bands in the 
energy space. 

According to the prescription of this approach, we 
pick out the one-electron band Hamiltonian Hb which 
describes the hybridization of band electrons with the 
mean- field level £/, which is defined as the energy differ- 
ence 

where Er-^ is the lowest (doubly degenerate) state in the 
multiplet f^(^iJ5/2), and then represent the Hamiltonian 
(|I|) in the form H = Hb + AH. The term AH describes 
the excitations above the ground state. Then the mean 
field spectrum in the actual energy interval can be de- 
scribed approximately by a Hamiltonian 



(2) 



= 1,2 k(T 



Here 



£i,2(k) = i(6k + Sf) T ^J(ek - Sff + \VikW, (3) 

where e"" ''V^(k) = V^k from Eq. (^. Hybridization 
between p- and /-states due to the matrix elements F (k) 
results in the wavefunctions 



(4) 



where the coefficients of the canonical transformation are 



Uk 



Ek 



(5) 



Since the band states and the /-states have the same 
symmetry at the points X7, R'j Fy and Ay a hy- 
bridization gap opens, and the lower band £1 is filled 
in the ground state of "divalent" SmBe because two 
electrons transferred from the /r^ level to the (initially 
empty) band ei fill it completely. Let us assume that the 
level £/ crosses the band e(k) in its lower part, so that 
the band ei(k) has mainly /-character. Then the square 
root in Eq. (^ can be expanded in F(k)/[ek — £/] for 
most of the Brillouin zone, and the Wannier states given 
by the operators ci^mtr = N^^^^ exp(jk • m)ci_k(T can 
be approximated by the following equations. 



Cl, 



7 ~ fma 

k (j)jviv 



^(m-j)cxp[zk- (m-j)] ^ 
ek-£/ 



(6) 



Here {i)NN are the nearest neighbours (boron sites) of 
a Sm ion in a given crystal cell and V{m — j) is the 
Fourier transform of V(k). The width T of the lower 
band is determined by a hopping integral which can be 
estimated as 



kj 



V{m — j)V^*(n — j) exp[zk • (m - n)] 



ek-£/ 



(7) 



(here j are the nearest neighbors of both m and n) . 

To study the real carriers, one should go beyond the 
mean-field Hamiltonian, i.e. take into account the term 
AH which includes, in particular, the interaction be- 
tween holes in the narrow band ei(k) and the electron- 
hole interaction between carriers in different bands. The 
former contribution stems from a strong Coulomb in- 
teraction between electrons in the Sm /-shell U = 
(/m/m|C/|/m/m)- This Hubbard interaction is reduced 
due to slight delocalization of Wannier functions 

U = (ci,mCl,m|C/|ci,mCl,m) = r'^U, 

where r = J S{e)u{e)de is the reduction (nephelauxetic) 
factor j22j. This factor depends on the density of states 
S'(e) in the unhybridized band ek (see Eq. nf), and the co- 
efficient u{e) is that which appears in Eq. ^) but written 
as a function of energy. In our case the condition T <^ U 
is assumed to be valid, and this means that no more than 
a single hole per site can be created in this narrow "Hub- 
bard" band. 

The upper band £2 (k) is formed mainly by p-electrons 
with admixture of an /-component in the states close to 
the bottom of this band. If the bottom of the conduction 



band is close to the center of the Brillouin zone |20|, the 
Bogolyubov's coefficients for the electron wavefunction 
described by the operator C2.w M) are given by 



1 - 



Vo 



Uk 



£f - Ek 



(8) 



The Wannier operators C2,m have the same form as §), 
except that the sign of the denominator in the second 
term is negative. The dispersion in this band can be 
approximated by the effective mass law 



£2{k)^Ao- 



ri^{k-kbf 



(9) 



where Aq is the gap in the two-band spectrum (^ . Here 
kb is the wavevector corresponding to the bottom of the 
empty band, and the effective mass can be estimated as 



m:/mo «[(£/- 6feJ/y(fc,)]2 = 



(10) 



The effective mass is noticeably heavier than the bare 
mass mo of the conduction electron (experimentally, 
m* ~ lOOmo ||]). We suppose that the hybridization 
gap is seen as the biggest gap Aopt in the experiments 
mentioned in the Introduction. 
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Now, taking the energy of the state with filled lower 
band as a reference point in our calculations of the excita- 
tion energy spectrum, we have the effective Hamiltonian 



mn.cr ma 
+ £2(fc)4,k^C2,ka + iJl2 



(11) 



kcr 



where tmn — ~^^iO^) ^^'^ £i,2(k) are given 

by Eq. The electron states in the hopping term 

for hole excitations are dressed in projection operators, 
ci,mcr = Cl_mo■?^l,m-cr• This projcction inserts the kine- 
matic restriction for the hole motion: the hole can be 
created at a given site only provided the electron with 
the opposite spin is still at the same site. Therefore, 
the holes are practically immobile. The last term H12 is 
responsible for the electron-hole interaction, 



same crystal point symmetry as an /-electron in the cen- 
tral cell . This state is constructed from electron-hole 
pairs ( p^ by means of an envelope function Fq(m — j). 
When q = 0, the exciton operator can be written as 

\^e.)=N-'/'J2 E Foi^-i)pUiJma\0) . (14) 
m (j)jV]V 

Here a basis of / and p-orbitals is chosen. As mentioned 
above, the /-electron has r7 symmetry, and its angular 
dependence is determined by the xyz cubic harmonics. 
The envelope amplitude Fq is a function of exciton energy 
Eex (see, e.g., Ref. |l^] where this function is calculated 
in the NN-approximation) . It also contains the phase 
factor which orders the phases of the p^, Py and pz boron 
orbitals in accordance with the required symmetry (cf. 
II), 



Hl2^ E E^('^l'''2,k3,k4)4 ki^Ci,k2a4,k3<T'C2,k4<T' 



E 

kik2k3k4 <j<j' 



(12) 

where W is derived from the Coulomb interaction [last 
term in the Hamiltonian (|l])]. 

The two-band Hamiltonian H contains hybridisation 
built into the electron and hole states, and the average oc- 
cupation numbers nuj = fif are formally less than one in 
the lower "Hubbard" band ei. However, it should be em- 
phasized that this deviation from integer value {hf — 1 
corresponds to the Sm^^ state) is not the true IV state. 
The one-electron picture given by the mean-field Hamil- 
tonian Hb (H) implies that interband transitions given by 
the operator 



Sn 



E4 



2,k+qcr'^l,ko- 



-iq m t 



2,m+jCT''l,m(T 



(13) 



form the fundamental branch of elementary excitations 
in this semiconductor. However, free electron-hole exci- 
tations alone cannot explain the unusual properties of IV 
SmBg. We believe that in the IV state an extra branch of 
charge transfer excitations exists. These are valence fluc- 
tuations, and are responsible for the unusual low-energy 
electron spectra and the numerous anomalies of the phys- 
ical properties of SmBg. According to the scenario sug- 
gested in |l^,^ and verified in [|ll|, the true IV ground 
state arises as a result of admixture of singlet excitonic 
states to the ground state of the same symmetry, and this 
admixture is non-negligible when the binding energy of 
the exciton is comparable with the band gap, i.e. when 
the system is close to the excitonic instability. To realize 
this scenario, one should construct the singlet exciton for 
the specific case of SmBg. 

According to the theory of IV rare-earth semiconduc- 
the singlet exciton is a bound state of a hole 



10 2, 



tors 

in the samarium /-shell and an electron spread over the 
p-states of the surrounding NN boron atoms with the 



Fo(m-j) = F(Se,)(-l) 



M„ 



(15) 



where the phase (— 1) j follows the arrangement shown 
in Fig. g. 

In fact, it will be shown below that the IV ground 
state resembles in some sense the Zhang-Rice (ZR) sin- 



glets 1 15 which are believed to be formed in Cu-0 planes 
of high-Tc materials. The Emery Hamiltonian |Q, which 
is the starting point for the description of hybridiza- 
tion between weakly interacting oxygen p-electrons and 
strongly correlated copper c?-holes is similar to the Ander- 
son Hamiltonian describing the hybridization of nearly 
free 6-electrons in the conduction band with strongly lo- 
calized electrons in samarium 4/-shells. The difference is 
that the ZR bound states of a local Cu spin and a hole 
distributed over p-orbitals of the surrounding O ions are 
formed as excitations in doped oxicuprates, while the IV 
singlets are formed in the ground state as bound state 
consisting of a hole in the /-shell of the Sm ion and an 
electron distributed over the p-orbitals of the surrounding 
B ions of the cation sublattice . The binding mechanisms 
are the antiferromagnctic dp-exchange in ZR case and the 
Coulomb /p-attraction in our case. 
Now, inserting ( p^ into (p^), we have 



-1/2 " 



, -Pmcr/mcr |0) 



(16) 



where is a localised state which is defined as follows. 
SmBg has a simple cubic lattice structure, with an Sm 
atom in the centre of each cell, and a Bg group on each 
corner. To simplify the treatment, we approximate the 
electronic wavefunctions for the boron clusters by the p- 
orbital directed along the diagonal of the cube (see Fig. 
|l|). These orbitals are responsible for the covalent bond- 
ing between the boron clusters and Sm sublattice |6j. 
We introduce the linear combinations 

pf^^-^(p.(j)+p,(j)+p.(j)) 



Pi^' = -^{Pxij) - Pyii) + Pz{j)), 
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etc. which represent p-type electronic wavefunctions cen- 
tred at j and oriented in the direction x+y+z. In analogy 
with the Zhang-Rice construction, we write down a lo- 
calised state for the Sm site m consisting of eight nearest 
neighbour p-orbitals oriented to have the same symmetry 
as a central /-orbital, as shown in Fig. |l|. 
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(17) 




FIG. 1. A 4/-orbital located at the centre surrounded by 
eight p-orbitals on the corners. The signs of the 

p-orbitals are chosen to have the same symmetry as the 
/-orbital. 

States on neighbouring sites are non-orthogonal be- 
cause of shared p-electrons. Therefore, it is useful to 
make a canonical transformation to the orthonormalized 
Wannier or Bloch operators d^, 

P^^Y. - = N-'^'- J2 /^y'^e-^'--, (18) 

m k 

where 



(19) 



^ X . 2 • 2 ^'z , .2 '^x 2 '^V • 2 '^z 

COS — sm — sm V sm — cos — sm — 

2 2 2 2 2 2 



. 2 ^x ■ 2 2 ^z 

sm — sm — cos — 
2 2 2 



(20) 



Then the exciton operator acquires the form 



mcT 



/m.|0) 



(21) 



Here the index i ~ 1,2... enumerates the coordination 
spheres, and the coefficients A; fall rapidly with increas- 
ing distance. 







A 


Ao 


(0,0,0) 


0.935 


Ai 


(1,0,0) 


-0.104 


A2 


(1,1,1) 


0.056 


A3 


(1,1,0) 


-0.052 



TABLE 1. The coefficients A(m — n) where m — n = 
ix + jy + kz. 

By construction, the states |0) and \^ex) belong to 
the set of eigenstates of the Hamiltonian H (|ll|) with 
nonlocal terms Hih = J^m^n^a^^nncl^^^ci^na and SH12 
(defined below) excluded. 

However, these nonlocal interactions are responsible 
for forming the IV state. In a simplest approximation we 
restrict ourselves by considering only the first term ^ Aq 
in the expansion ( pi| ) and take the nonlocal interaction 
in the form 

6H,2 - - ^)fLfnafLdma + H.C. (22) 

mncT 

This component of electrostatic interaction violates the 
point crystalline symmetry and induces fp hybridization 
at the site m in presence of a hole in the neighbouring 
cell n. The operator 6H12 has a property 



= T4^(m-n)|0„ 



where On stands for the cell n in the ground state con- 
figuration. Thus this interaction intermixes the states 
|0) = Hn I On) and l^ex)- The mixing constant is u; = 
{0\SHi2\'^ex) — zF{Eex)^oW whcrc z is the coordination 
number for Sm sublattice and W is the NN interaction 
matrix element. After diagonalization the local neutral 
states in each cell m are represented by the linear com- 
binations 



\0m) =cos6l|Om)+sin6l|Vm), 
iV'm) = - COS 9\^m) + sin 9\0m 



(23) 



where tan 2^ « 21x1 /E^x^ and the valence is determined 
by the value of sin^ 9. As was mentioned, the local states 
( p3| ) are in some sense the electron-hole analogs of two- 
hole Zhang-Rice singlets and triplets (l^j2^ . Finally, the 
ground state of IV semiconductor is 



I*, 



(™)\ 



(24) 



and the low-lying local excitations are described by the 
vector 
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N- 



|V^m)(Om|* 



(25) 



The valence and the energy scale of these local ex- 
citations (valence fluctuations) are characterised by the 
degree of admixture of the excitonic state, in which an 
electron is promoted to a loosely bound "molecular orbit" 
Pm- If the exciton energy E^x is small enough (the ex- 
citon binding energy is comparable with the gap width) , 
the value of sin^ 9 is close to 1/2. Having in mind also 
the mean- field part of hybridization given by Eq. (^) , we 
conclude that in this case the valence 



Hf 



(26) 



can exceed 2.5 both in SmS and SmBg. Thus, the large 
deviation of valence from the integer value and the soft- 
ness of local valence fiuctuations are, apparently, the cor- 
related phenomena. Further discussion of internal con- 
sistency of the model can be found in the last section. 



III. TRAPPING OF CONDUCTION ELECTRONS 
AND HOPPING CONDUCTIVITY 

The samples of good quality which were studied in the 
experiments of the past decade mentioned in the Intro- 
duction are n-type semiconductors at low temperature, 
unlike the p-type samples of the first generation |2g] . The 
electron concentration in these samples is estimated at 
Ue ~ 10^^ cm"'^ at ambient pressure and liquid helium 
temperature |^,^. In order to interpret the transport 
properties of these samples one should determine the 
spectrum of electrons at the bottom of the conduction 
band in the presence of soft valence fluctuations. In our 
model the local valence fluctuations arise as transitions 
between the states |0m) and \'tpm) given by Eq. (p3|). It 
is known that these valence fiuctuations are the source 
of strong anomalies in the vibration spectra because the 
characteristic time r^f of valence fluctuations is close to 



phonon times Tph 



OJ, 



ph 



10 s. Therefore, one can 



expect that these "slow" excitations could dress the car- 
riers and form an electron-polaron cloud similar to the 
phonon cloud which results in polaron self-trapping in 
dielectric crystals (see, e.g., [p7|). 

To describe electron self-trapping we start with a 
Hamiltonian including the interaction between the con- 
duction electron and the local valence fluctuations in a 
single lattice site n = 0, 



i/e = ^e(k)c[ck + r!o^^^ 

k 

[l^./(kl,k2)4^Ck,At+i/., 



(27) 



ki,k2 



This Hamiltonian stems from our basic Hamiltonian (11). 
It includes the electrons in the upper conduction band 
(the spin summation and band index s = 2 are omitted) . 
Only the term which corresponds to momentary, local 
redistributions of charge and valency at the site of the 



excitation due to the interaction with the charge carrier 
is retained in Hi2. The valence fluctuations with the en- 
ergy rio are described by the local operator = |V')(0|. 
This excitation is in fact a charge transfer between the 
periphery and the centre of the cell. The "breathing" 
mode r^o describes the polarisation of Sm ions in the 
cation sublattice that accompanies the propagation of 
an excess conduction electron (predominantly, over the 
samarium sublattice, see below). 

The matrix element W^/(ki,k2) is given by the inte- 
gral 

W^„/(ki,k2) = 

dridr2V5*^(ri)?/;o(ri)VF(ri,r2)V'ki(i"2)-0k2(r2) 

where '0k(r) are the Bloch functions for the conduction 
electrons, "iplxi^) and '0o(r) are the wavefunctions of the 
states created by the operators (^). Their product is 

r,x{r)Mr) ^ D{v) 

= i (sin20[pp(r) - pf{Y)] - COS29 fxyz{v)dxy^{Y)) , (28) 

where [r) — \fxyz{f)\'^ is the f-clcctron density in 
the centre of the cell (see Fig. 1) and /9^(r) = 
F'^ {Eex)^Q\dxyz{r)\^ is the p-electron density in a cen- 
tral cell given by the first term in the r.h.s. of Eq. (|2l]). 
The conduction electron density operator /5ki,k2(ri) = 
■0kj (r2)^k2 (r2) is determined by the Bloch functions (|^). 
Near the bottom of conduction band this operator can 
be approximately presented as 



Pki,k2(ri) ~ {r) + CcpI,m.,{yi) 
by means of Eq. (||) . Here 



Pi 



(S)nn 



(29) 



(30) 



(31) 



and Efb = ~ £feb- Then the coupling constant in the 
Hamiltonian (|27| ) acquires the form 



W„/(ki,k2) » wo + u;i/3ki-k2 (32) 



where 



Wo = J dridr2D{ri)W{ri,r2)pHr2), 
wi/3ki-k2 ^ J dridr2D{ri)W{ri,r2)pl^ i,^ir2 + 3). 



The structure factor /3k is given by 



> cos — cos — cos — 
2 2 2 



(33) 
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It is seen that the first term wq dominates in the electron- 
exciton interaction, therefore we begin with a calculation 
of the electron self-energy induced by this term within the 
lowest order of perturbation theory. The theory which al- 
lows inclusion of nonlocal corrections due to small term 
^ wi in Eq. (|3^ ) is described in Appendix. We calculate 
the electron Green's function in Matsubara representa- 
tion 



rl3 

Gkk'K)=/ {TrCkiT)cl)e'^"^dT 
JQ 



(34) 



[ujn = (27^ + l)7rT] interacting with a localized valence 
fluctuation mode Qq described by the propagator 



2^0 



(35) 



[uj^ = 2'KmT]. The Dyson equation for the electron 
Green's function is 



Gk,k'(w„) = G£(t^) 



^k,k' + ^ Sk,k2('^ri)Gk2,k'(ti'n) 
k2 



(36) 



In the lowest order of perturbation theory the self-energy 
5jki,k'(wn) after carrying out an analytic continuation 
into the region of real frequencies acquires the form (see 
Fig. 2) 

S^,k2H = E^-/('^'l^i)^-/('^i'1^2)7'ki(w,r!o) (37) 

ki 

where 



J-oo 2Tr 



e + e(ki) " u) — iS 
2T 



tanh -—^ coth — 



2T 



(38) 




FIG. 2. The self-energy diagram corresponding to Eq. 
(37). The straight line represents a conduction electron and 
the wavy line represents a soft valence fluctuation. 

Approximating the nonlocal potential VFu/(k, ki) by 
its local part wq, we come immediately to the following 
secular equation 



1 — WQm(ti>)m(a; — ilo) = 



(39) 



where miuj) = X]kGk('^)- 



In close analogy with the polaron self-trapping effect 
p7| , one can expect that the attractive polarisation po- 
tential with effective coupling constant 
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results in appearance of bound electron-polaron states. 
As usual in 3D problems, the threshold value of the 
attractive potential determines the onset of the bound 
state. To find the constraints on the solutions with nega- 
tive energy (corresponding to bound states) we examine 
the function m{uj) defined in (|3S|), 



m{uj) = 



dkx dky dk^ 



1 



(27r)3 w-e(k)' 



(40) 



At a; = the integral may be evaluated numerically (see 
Appendix) . Fig. 3 demonstrates the graphical solution 
of Eq. (|3^) for a reasonable set of model parameters. 
The value of JIq = 5.5meV correlates the experimentally 
observed peak in optical reflectivity spectra [g|j7| , and the 
binding energy of localized uji ~ 3.5meV is in good agree- 
ment with the activation energy registered in multiple 
optical and transport measurements in the temperature 
interval 6 - 14 K (see |^ and references therein). 

We have demonstrated the existence of a self-trapped 
state in the simplest approximation, whereby the local 
polarisation mode is taken into account in second or- 
der perturbation theory. However, the values of the pa- 
rameters necessary to achieve reasonable agreement with 
experiment give for the dimensionless coupling constant 
a = wo/^lo the value of a « 4.5. This means that in 
fact we are in a strong coupling limit, and a more refined 
treatment is necessary. The generalisation of the theory 
to the strong coupling limit can be done in close anal- 
ogy with the theory of small polaron. One can make the 

where 



canonical transformation Ck = e '^CkC'^ 



kq 



W^./(k,q) 



C^Ck+qA^ 



4+qCkA_q) (41) 



(see, e.g., [^). We have seen above that the local 
term wq is dominant in electron-exciton coupling con- 
stant (32). Then neglecting the contribution of the "tail" 
wi we come to purely local interaction at a given site 
m = 0, namely, SH12 = woCqCo{Aq + Aq). Eliminating 
this interaction by means of a canonical transformation 
(^), one comes to the effective Hamiltonian 

He = e-^^H^e^ 

m m, 115^0 

- e^oiclco - J2 {Tmocic^e-«-^°) +h.c.) . (42) 

m=iO 

Now the main part of electron-exciton interaction is 
taken into account exactly, and it is contained in the po- 
laron shift. This shift induces the local scattering at the 
site which is given by the operator Epd k2 ''ki''k2J 
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where Spoi = aw^. The local scattering can be inserted 
in the self-energy of the conduction electron in the same 
manner as it was done in second order of perturbation 
theory in Eq. (36). However, now the calculation is ex- 
act, and instead of Eq. (B9h one has 



= l-f Epoi ^[^- e(k)] ^ = 1-1- epoim{uj) (43) 

k 

which can be solved by using the same approximation as 
in Eq. (|3^). In any case the electron can be trapped 
by local valence fluctuations provided the polaron shift 
£poi — awQ given by the bound solution of Eq. ex- 
ceeds the characteristic kinetic energy of the electrons 
near the bottom of the conduction band. Therefore 
the effective mass enhancement due to ^/-hybridization 
favours the formation of a self-trapped state. 



2 - 



- 2 



3 
3 




3 



Cv? C 

if 



cj(tti eV) 

FIG. 3. Graphical solutions to Eqs. (39) (solid line) 
and (43) (dotted line) using eo = 5meV, m* = lOOmo, 
f^o — 5.5meV, wq = 25 meV and Spoi = 21meV. 

In Fig. 1^ we show a bound state solution with en- 
ergy 3.5 meV in agreement with the activation energy ob- 
served in transport measurements and absorption edges 
in optical experiments. In both cases the coupling con- 
stant is treated as a free parameter. In the second 
order perturbation calculation wq is determined to be 25 
meV, which as we have discussed above, indicates that 
our solution lies in the limit of strong coupling. In the 

non-perturbative calculation we find that epoi = ^ = 21 
meV, which implies that wq = 11 meV. Because we are 
in the strong coupling limit we do expect a deviation 
between the exact calculation and the second order per- 
turbation result. However, our results are close enough 
to suggest that the physical content of the perturbative 
calculation is correct, that is, that an electron may bind 
to a single local valence fluctuation and the result is a lo- 
calised electron-valence fluctuation complex. Moreover, 
we also observe that our results are not particularly sen- 
sitive to the input parameters, and are therefore not the 
result of any special "flne tuning" . All of this is a con- 
sequence of the fact that the relevant energies are all 
roughly the same: VIq = 5.5 meV, h^/a^m* = 5 meV 
and Aact = 3.5 meV. 



The last term in the effective Hamiltonian (42) is re- 
sponsible for the exciton-assisted hopping in an electron- 
polaron band. The electron self-trapping can occur 



around any arbitrary lattice site, so the trapped parti- 
cle can move from a given cell to the neighboring cells 
by the mechanism which resembles the polaron propa- 
gation |2^,|2^. Therefore, the electron-polaron drift in 
electric fleld is responsible for the conductivity at T ^ 
in n-type SmBe. We leave systematic treatment of elec- 
tron transport in IV semiconductors for future publica- 
tions, and conclude this article by qualitative discussion 
of anomalous low-T behavior of SmBg mentioned in the 
Introduction. 



IV. CONCLUDING REMARKS 

As a result of above analysis, one can understand the 
mechanism of reduction from the eV energy scale char- 
acteristic of initial Hamiltonian (|l|) down to meV scale 
of low-lying charge excitations in the n-type material. 
First, at the mean-field level of approximation (eqs. ^ to 
|lO|) the gap Aq ~ 10 — 20 meV opens, and the states near 
the bottom of conduction band become heavy. This ini- 
tial reduction of the energy scale is characterized by the 
parameter Q (|30|), and the position of the /-level relative 
to the bottom/top of the bare conduction band ek can be 
extracted from the ratio of effective masses 



where eft = £kt - 



1/2 



1/2 



(44) 



= tfet — £/, f-kt is the top of the band e^, and 
the parameter = (Vb/£/t)^ is determined similarly to 
Cc- We neglect in these crude estimates the anisotropy of 
hybridization integral T^(k). The hole effective mass m* 
is evaluated as ^ 500— IOOOtoo [§j5|, so we take r ~ 1/3, 
which is consistent with our assumption that the mean- 
fleld hybridization gives only part of the actual value of 
the valence According to our previous calculations 
p2t , the value of r = 1/3 corresponds io nj 0.75, and 
the experimentally observed value of n^, « 2.55 can be 
reached only by means of excitonic mechanism with the 
exciton mixing parameter svnO ~ 0.55 (see eq. ^6|). This 
value, in turn, agrees with our assumption of a soft ex- 
citonic mode. Then, we find from (|l^) Qc — 0.01, and 
taking for the bandwidth Sfb + £ft the value of 4 eV in 
accordance with the band calculations for the related sys- 
tem Q , we estimate the hybridization coupling constant 
as Vq « 0.1 eV which is a reasonable value for the rare- 
earth materials. To make these estimates self-consistent, 
the hybridization gap should be found. The gap is de- 
fined via the above parameters as Aq « Cc£/b + Cv£ft- 
Inserting the values of corresponding parameters, we find 
Ao ~ 13 meV which is in a reasonable agreement with 
the experimental data. 

Now, turning to the n-doped materials, we deal with 
heavy electrons, which interact with the valence fluctua- 
tions. The kinetic energy of these electrons is estimated 
as ^ Vo/£fb ^ 10~^ eV , and the energy scale of va- 
lence fluctuations is, by their origin, limited from above 
by the energy gap of the same 10"^ eV width. The po- 
larization coupling constant Wyf (32) as well as the "su- 



pcrhybridization" matrix elements W in the Hamiltonian 
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( P4 ) should be at least an order of magnitude less than 
the mean-field hybridization V, so the value of wq = 25 
meV, which is used in our numerical solution of Fig. 3 
looks realistic. Eventually, solving the equation (^3|) for 
the polaron shift, we descend one more step along the 
energy scale and find ourselves in a meV region. The 
propagation of self-trapped electron-polaron can in prin- 
ciple be characterized by even lesser energies of the order 
of 10-1 meV. 

Thus, we have found that the heavy electrons near the 
bottom of conduction band can propagate in a lattice 
only in a polarization cloud of valence fluctuations. It is 
worth mentioning that this conclusion correlates partly 
with a recent proposition of Kasuya p9[ | in spite of the 
fact that his model of the ground state of a mixed va- 
lence semiconductor disagrees with our picture in several 
respects. He chooses the trial wave function for a local 
singlet (which he calls Kondo singlet) of the form 



|0m) = a/m|Om) + adm^rfml'^n 



(45) 



where ^j^^ = dl^Cdm is a charge transfer operator, which 
creates an electron in a 5d shell of Sm ion and a hole of d- 
symmetry spread over surrounding boron sites similarly 
to our f-like orbital entering the state |V'm) (pl|)- The 
physical reason of such a choice is conviction that the 5d 
electrons also constitute a strongly correlated subsystem 
and have a trend to form a bonding Kondo-like state a la 
/ states in conventional Kondo lattices. However, looking 
at the real band structure of the rare-earth hexaborides, 
we see that the admixture of d states to the band ek is 
really small: the center of gravity of d states in divalent 
hexaborides is at the energy of ~ 7—8 eV above the bot- 
tom of conduction band (see, e.g., Fig. 5 in Ref. [pTf). 
Therefore, 5d levels are practically unoccupied, and one 
cannot expect any kind of effective screening in a d chan- 
nel. Nevertheless, our model includes effective screening 
of excess electron in conduction band which resembles 
Kasuya's mechanism: the canonically transformed oper- 
ators Co in a weak coupling limit can be presented as 
Co ~ co(l -I- Q!^q), and the operator returns part of 
the charge density back to the Sm site from the periphery 
of the unit cell, i.e. plays the same role as the operator 
yljjQ in the trial function (^5|). The screening (or po- 
laron dressing) in our model is due to the same valence 
fluctuations as all other physical effects, so the totality of 
experimental data is explained in a self-consistent scheme 
without appealing to any additional hypotheses. 

As a result, our model is free from undesirable features 
of "Kondo insulator" approach which are not confirmed 
by the experiment. 

(i) We do not appeal to the Kondo mechanism of forming 
the ground state singlet, and the independence of the ac- 
tivation gap on external magnetic field [Q, which rules 
out the Kondo insulator mechanism, agrees well with the 
expectations of our model. 

(ii) The ground state of our Hamiltonian is absolutely ho- 
mogeneous, and this statement agrees fairly with avail- 
able experimental observations, whereas the trial func- 
tion (Bq) implies charge modulation in a form of Wigner 



crystal or Wigner liquid 

(iii) The model of ref. |28fP gives a single gap which pre- 
determines electronic, optical and magnetic properties of 
the material, and it is unclear whether it is compatible 
with a real experimental situation which definitely ev- 
idences several energy scales for electronic and optical 
characteristics of SmBg |^,||- 

(iv) Unlike the bare ground state |0), the wavefunction 
( p^ ) is not fully symmetrical, so one can expect a sort of 
ferroelectric ordering as T ^ ||3l|,^ . 

Turning to the problem of low-temperature transport, 
we conclude that as T ^ the electron propagation is 
valence fiuctuation assisted motion of in extremely nar- 
row "polaron band" separated by the gap Aact from con- 
duction band continuum. The residual resistivity po is 
inversely proportional to the hopping rate between two 
neighboring crystal cells. 



Po 



Sr, 



T $ 



(46) 



where $rnn — ("0111 IV'n) is a function describing the over- 
lap of the valence fluctuation "clouds" centered around 
the sites m and n. As T — > the process is elastic, so the 
conductivity is temperature independent. In SmBg this 
type of conductivity is observed for T < 3K ^,||,D. At 
higher T the hopping is, apparently, assisted by phonon 
and exciton emission/absorption. Of course, the temper- 
ature dependence of this hopping can differ from that for 
usual variable range hopping in impurity bands. This 
transient regime can be seen in the range 3 K < T <6 K, 
although the experimental data on the temperature de- 
pendence are still ambiguous [^|j5|J|] . In the temperature 
interval 6 K < T <14 K thermally activated resistivity 
in SmBg with an activation energy of Aact ~ 3.5meV is 
observed At higher temperatures the electrons 

dissociate from their local valence fluctuation clouds and 
as a result of this detrapping flnd themselves at the bot- 
tom of the conduction band. At these temperatures the 
activation of valence electrons also gives a significant con- 
tribution to the electron conductivity. 

The same three regimes (III, II and I, respectively, 
in terms of Ref. |^) with additional pronounced maxi- 
mum around 5 K are observed in the temperature be- 
havior of the Hall constant Rh{T) This maximum 
can be explained, at least qualitatively, within a sim- 
ple phenomenological picture of two groups of carriers 
with high and low electron mobilities, /^f, — cRbUb and 
/i/i = cRh<Jh, where Ri, h and ab,h are the Hall constants 
and conductivities of light (5) and heavy (h) carriers, re- 
spectively. In the case of hopping in impurity bands of 
doped semiconductors two contributions, Uh (hop- 
ping) and (Tft (band) in the electron conductivity result 
in the following equation for the Hall constant, 



Rh 



Rh<jl 



(O'b + CThY 



(47) 



Then the maximum in Rh (T) corresponds to a crossover 
from hopping motion at low T to band motion at high 
T, provided fib^ fJ-h- 
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In our case the phenomenological background of this 
equation still exists, but the microscopic origin of all tem- 
perature dependences should be revised because of the 
essentially many-particle nature of heavy carriers. 

First of all, it is clear that the standard estimates of the 
number of "scattering centers" obtained from the value 
of the residual resistivity po are simply inapplicable in 
our case, since we deal with hopping of many-particle 
electron-exciton complexes rather than with the motion 
of extended band electrons. So, there is no room for the 
unitarity limit arguments in these estimates, and the "su- 
perunitarity scattering" reflects the many-particle nature 
of current carriers in semiconductors with fluctuating va- 
lence P^L Thus, the paradox of the number of scatterers 
per site|^ is removed, and we return to the usual situa- 
tion with one scatterer per unit cell. 

The next question is the enormously strong pressure 
dependence of the residual resistivity and the Hall coef- 
ficient jij. Again, according to our model, one cannot 
directly apply the notions of charge transport in a band 
of extended states for an estimation of the carrier con- 
centration. We think that the key to the extraordinary 
sensitivity of the residual resistivity po and Rh (T — > 0) 
to external pressure is the increase of Sm valence with 
growing pressure P. The eventual source of this increase 
is growing /p-hybridization. Because of increasing va- 
lence function of pressure the matrix clement 
Wyf (^) decreases. As a result the polaron shift £poi 
as well as the binding energy of the self-trapped electron 
decreases too. The excitonic overlap function <&inn f^lso 
grows due to lattice contraction. Since both the hopping 
integral Tmn and the function <I>mn depend exponentially 
on the intcrsite distance, one can expect a very sharp de- 
pendence of the hopping rate on external pressure. The 
radius of the localized state with binding energy 3.5 meV 
and effective mass ~ lOOmo is estimated to be 2-4 X 
and an increase of this radius by an order of magnitude 
would be enough to get a 10'* growth of the hopping 
rate. Such an increase is achievable under a pressure of 
50 kbar. At higher pressures the trap becomes too shal- 
low (the valence too close to 3) to catch the electron, the 
system transforms into conventional n-doped degenerate 
semiconductor, and eventually it becomes a metal with 
trivalent Sm ions in the cation sublattice. 

Finally, the question of the variation of electron con- 
centration rie with increasing P and T also demands spe- 
cial consideration. Usually information about is ex- 
tracted from the value of Rh under the assumption of 
single band conductivity when Rh « i?f, — (tIeCc)^^. 
According to the experimental data cited above, Rh is 
nearly constant in the temperature interval I (below 3 
K), but falls drastically with applied pressure: i?//(45 
kbar)/i?fl'(l bar)~ 10"''. However, in our case of two- 
component systems with ^ cr;, the Hall constant takes 
the form 

Uh ~ 9 . 

It is seen immediately that the pressure dependence of 
this function is determined mainly by the denominator, 



and the decrease of Rh with pressure is due to the ex- 
ponentially growing factor S^n in Eq. (||) for fffi rather 
than to increasing Ue- 

Generally speaking, Rh derived from Eq. ( ^7| ) also 
cannot be used to determine electron concentration nei- 
ther in region I nor in in region II. The low T limit of 
rteo ~ lO'^^cm"'^ which is obtained under the assumption 
that Rh ~ ^6 can only be considered as a lower bound. 
In the general case of two mechanisms of charge trans- 
port the Hall constant is Rh = {ni,ec)~^ p,H / fJ'Di where 
fJ'H.D are the Hall and drift mobilities respectively. If 
tJ^ah^D S> 1 at low T the real electron concentration n^, 
may substantially exceed the value of n^Q . Referring to 
the exponential dependence of <t{T) and Rh{T) in re- 
gion II we should note that there are at least two causes 
of such a dependence. First there is the thermal acti- 
vation of electrons from the polaronic traps to the band 
continuum states, and second is the temperature depen- 
dence of the hopping rate S'mn due to the contribution 
of thermally activated exciton assisted terms which ap- 
pear instead of the overlap integral <i>inn- As a result, 
both <Jfi{T) and (Tfc(T) give an exponential contribution 
with the activation energy of Aact ^3.5 meV to the elec- 
tron conductivity and the Hall constant ( ^6| ) . Again, it is 
impossible to calculate the variation of the electron con- 
centration ne{T) directly from this equation. We leave 
the detailed evaluation of various transport coefficients 
for a forthcoming paper, but to conclude this qualita- 
tive discussion we would like to emphasize that accord- 
ing to the theory proposed in this paper the total carrier 
concentration in conduction band is ~ lO'^^ cm"'^ at 
least while the interband electron activation is negligible 
(T <C Aq), and there is no room for dramatic pressure or 
temperature variation of in the mixed-valence phase 
of n-SmBg. The real semiconductor to metal transition 
occurs only at P > Pc when the Sm valence changes to 
the integer value of -1-3. 

The authors arc grateful to M. Aronson, D.I. Khomskii 
and N.E. Sluchanko for useful discussions. This work was 
supported by Israeli Academy of Sciences and Human- 
ities Centre for "Strongly Correlated Interacting Elec- 
trons in Restricted Geometries". S.C. acknowledges the 
support of the Feinberg School of the Weizmann Insti- 
tute. 



APPENDIX 



We describe below the self-trapped state of the elec- 
tron captured by the nonlocal potential (|3^). For the 
sake of simplicity we consider the case of conduction 
band with a minimum at the T point of the Brillouin 
zone. To solve the Dyson equation ( ^6|) we offer the 
procedure of factorization of the vertex matrix element 
W^/(ki,k2). The latter can be represented in a form 
I^.^(ki,k2)^ELo7S,7£„ where 
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7k 



y/8(wi COS ^ COS ^ COS ^ 

y/8(wi sin % COS ^ cos ^ 



- 2 

V8Ciyicos^sin%cos ''^ 



"2 '^J^-ii 2^ 

VSCtyTcos^^cos " 



2 

\/8^^t^j"sin ^ sin ^ cos ^ 
VSCwi cos ^ sin % sin ^ 
V8Cwi sin ^ cos ^ sin ^ 
V V8Ci^sin^sin^sin% 



(48) 



By writing G£ = X^k' 7£'Gk,k' 
for G of the form 



we obtain a solution 



(49) 



where 



M"" = -J''" - ^ ^ G°^ (a-)7'k. (^, f^o)7£,7£,7£,7L 



Here the index a = corresponds to the fully symmetric 
solution described approximately by Eq. (P9| ) , and a = 1 
stands for the nonlocal contribution of the "tail" wi with 
the same Ai symmetry. For the indices a,b = 2 ... 8 the 
matrix M is actually diagonal since terms in the inte- 
grand must be even in k\x^ ^ly? ^iz and k2x^ ^2^? ^22 and 
therefore a = c = b. For these components M is given by 



M"" = dab [1 - {10)771" {lo ~ no)] 



a, b 



(51) 



For the components a,b = 0,1 M reduces to a two by 
two matrix 



/ 1 — 7n^ (Lj)7n'^ (uj — Q,o) —7nP(Lj)m(uj — Hq) \ 
— TO(w)m(ci; — r^o) —171(^1)7)7,^ {lu — Hq) 



M 



—7fi{oj)mP {(jj — Q.o) 1 — m(w)m(a; — flo) 
\ —7n}{u))7h(u) — rio) —7n^{u!)m^{uj — flo)/ 



In Eqs. (51) and ( |52| ) m{u}) is given by 

m'^(c.) = ^G°(c.)(7£)2 a = 0...8 

k 

m(c.) = ^G£(a;)7£7^. 



(52) 

(53) 
(54) 



= DetM 

= [1 - m"(w)m°(tj - rio) - 277i^{uj)m\u; - Qq)] 



Jl [1 - m''(w)m''(w - ilo)] 



(55) 



a=2 



m 



, , /" dkj; dky dk^ 

C(l ± cosfc2:)(l ± cosA:y)(l ± coskz) 
Lo - e(k) 



(56) 



Here the choice of three -|-signs corresponds to the Ai 
state, two -l-signs corresponds to the Ti triplet, two 
—signs corresponds to the T2 triplet and three —signs 
corresponds to the singlet A2. To study the analytical 
properties for small w < 0, we use the fact that the dom- 
inant contribution to the integral comes from small values 
of k. Near the bottom of conduction band e2(k) the hy- 
bridization is strong, and the band is nearly flat [see Eq. 

)], so expanding the dispersion around the minimum 
2 1,2 , 



(50) at the bottom, e{k) « h k"^ /2m* in accordance with Eq, 
(^), we refer to the heavy effective mass of m* « IOOtoo, 
observed experimentally 1^,^. Then we find for small 
w < 0, 





0.39 



1.31 



-0.35 + 0.28— 

-0.20-1-0.05— 

\uj\ 

-0.14-1-0.02— 



where eo = It' /a 



(57) 

(58) 

(59) 
(60) 
(61) 
(HI) for 



' « 5meV. Solutions of Eq. 
the values of the parameters given in Fig. 3 yield no 
bound states for Ti, T2 and A2, so we conclude that the 
approximation W^j ~ for the potential ( |3^ ) is suf- 
ficient for the description of the bound electron-exciton 
states. 



Thus we find that M is reducible to the different rep- 
resentations of the group of symmetry operations on 
an octahedron Oh'- a = 0, 1 and a = 8 stand for the 
one-dimensional representations Ai and A2 respectively, 
while a = 2, 3, 4, and a = 5, 6, 7, are the indices of triplet 
states Ti and T2. 

Hence, the secular equation which generalizes Eq. ( |39| ) 
has the form (for ^ <C 1) 
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